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Abstract 

In the Standard Model of particle physics, the quarks and anti-quarks have fractional charge equal 
to ±1/3 or ±2/3 of the electron's charge. There has been a large number of experiments searching 
for fractional charge, isolatable, elementary particles using a variety of methods, including e + e~ 
collisions using dE/dx ionization energy loss measurements, but no evidence has been found to 
confirm existence of free fractional charge particles, which leads to the quark confinement theory. 
In this paper, a proposal to search for this kind particles is presented, which is based on the 
conservation law of four-momentum. Thanks to the CLEOc and BESIII detectors' large coverage, 
good particle identification, precision measurements of tracks' momenta and their large recorded 
data samples, these features make the scheme feasible in practice. The advantage of the scheme is 
independent of any theoretical models and sensitive for a small fraction of the quarks transitioning 
to the unconfinement phase from the confinement phase. 
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In the Standard Model (SM) of particle physics, the mesons comprise two of quarks and 
anti-quarks and the baryons comprise three of quarks and anti-quarks, while the quarks and 
anti-quarks carry fractional charge ±l/3e or ±2/3e, e being the magnitude of the electron 
charge' 1 '. The successes of SM in the fields of the spectroscopy and the high energy inter- 
action of the hadrons have motivated a lot of experimental searches for fractional charge, 
isolatable, lepton-like particles. There has been a large number of experiments searching 
for free fractional charge particles using a variety of methods and no evidence up to date 
has been found to confirm their existence. There are three kinds of experiments used to 
search for fractional charge particles, which are based on the collision products of accelera- 
tor beams' 3-6 ', cosmic rays' 7 ', and bulk matters' 8 ' respectively. No observation of fractional 
charge particles results in the confinement concept that the quarks are permanently confined 
within the hadrons. The questions are whether there is a small fraction of quarks transi- 
tioning to the unconfinement phase from the confinement phase and whether there are new 
sensitive methods to search for them. 

It is noted that dE/dx ionization energy loss measurements are used to identify charged 
particle species for e + e - ' 3 ' 4 ' (or proton anti-proton' 5 ' or heavy ions' 6 ') collisions in accelerator 
based experiments. The energy loss fluctuation is described by Landau's theory' 9 '. One 
wonders if the theory is still valid for quarks, which take part in both the electroweak 
interaction and the strong interaction and is sensitive enough to distinguish the quarks from 
the five stable particles such as the electron, muon ,pion, kaon and proton if there is a small 
fraction of quarks which are unconfined. However it is no doubt that the conservation law 
of momentum and energy (four-momentum) applies to all kinds of interactions, which can 
be used to search for fractional charge particles in combination with the Lorentz force law 
for low energy accelerator experiments. The prerequisite for this kind of search is that a 
detector is required to possess large coverage, good particle identification, precise momentum 
resolution of charged tracks and energy resolution of neutral tracks and collect large data 
samples. For instance, recently upgraded CLEOc' 10 ' and BESIIl' 11 ' detectors can search 
for fractional charge particles based on the laws. The center masses energy of the e + e~ 
collisions in CLEOc and BESIIl cover the energy range of the charmonium production. 
Therefore fractional charge particles in the mass range of to 2 GeV can be conceived to 
appear in the products of the collisions if they are physical reality. 

The momentum measurement of a track with charge qe produced in an e + e~ event is 
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based its bending radius R in the magnetic field B along z direction in the tracking chamber 
according to the Lorentz force law, 

P xy = qeB/R, 

where P xy is the transverse component of the track momentum. The longitudinal component 
is measured by considering the polar angle 6 of the track with respect to the e + beam 
direction, 

P z = P xy ctg(9). 

In other words, the radius is a measured value while the momentum is a derived value and 
dependent on the assumed charge q. In the CELOc and BESIII experiments' 10 ' 11 ' as well 
as others' 3 " 6 ' q is a priori set to be 1 no matter how large charge it has. The momentum 
P n measured in this way, we call it the nominal momentum, for the track equals to its real 
momentum P r if q is 1. If q is not 1, its nominal momentum changes by a factor of l/q 

— * — * 

compared with its real momentum, i.e. P n = l/qP T . 

Hereafter q is assumed to be less than 1 for simplicity. For events with at least one particle 
of charge larger than 1, the similar analysis can be performed. Assuming there is a pair of 
quarks / and / with opposite charge qe in the event such as ty + ti^uu, K + K~dd and ppss, 
which is required to have even number of charged tracks. The total nominal momentum and 
energy are derived as 

P n (tot) = P n (/) + P n (J) + others 

= 1/qPrU) + V<?£(7) + others , 
£ n (tot) = sJPAf) +mj + yJp^Q) + m) + others 

= y / (l/gP r ) 2 (/)+mJ+ ^{l/qP r f(J) + m) + others. 

— * 

If P n (tot) 7^ 0, there are two cases. One is that some tracks in the event are undetected for 
they are neutrinos or go beyond the coverage of the detector. Another is that q is not 1. 
For the latter case, q can be adjusted so that the real total momentum 

P r (tot) = qPn(f) + qPn(f) + others = 0. 

Then the mass mj of the quark (anti-quark) can be calculated by requiring the real total 
energy 

£ r (tot) = ^(qP n ) 2 (f)+m}+ ^{qPn) 2 (f)+m} + others = ^s, 
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there y/s is the mass center energy of the e + e~ collisions. In practical measurement and 
calculation, one does not know which pair of particles with opposite charge have fractional 
charge. The momenta and mass of any one pair of particles with opposite charge can be 
adjusted so that the event satisfies the conservation law of momentum and energy. It should 
be noted that P r (tot) is a vector with three components, one cpmponent can be used to 
derive q, the two others can be used to suppress events from the first case because they 
can not satisfy the conservation law of three components simultaneously while events with 
fractional charge particles can. The procedure is repeated for all events in the data sample 
and the parameter set of (q, m/) is obtained. Then the parameter set is plotted into a two 
dimension distribution. For the first case, (q, mf) are a set of random numbers, but for the 
latter case, they should concentrate at one or more points if there are one or more fractional 
charge particles in the data sample. If fractional charge particles carry continuum masses, 
they should concentrate along one or more lines. 

To suppress the contribution to (q, m,f) from the first case further, the velocity v m of 
each of the pair of particles, which is measured by the detector, can be compared with the 
derived velocity v r = qP n / \J m2 f + (<?-Pn) 2 - They should be consistent, i.e. v T — v m — if 
it is a real fractional charge particle. If the parameter set of (q, v T — v m ) is plotted into 
two dimension distribution, the points corresponding to fractional charge particles should 
concentrate along q axis. 

To search for events containing a pair of quarks fl and f2 with opposite charge qe and 
different masses m/i and mj 2 such ds, the events are required to have even number of 
charged tracks. The total momentum, energy and one of the velocities 

P r (tot) = qPn(fl) + qPn(f2) + others = 0, 
£ r (tot) = y / (gPn) 2 (A)+™/i + ^ ( q P n f(J2) + m 2 j2 + others = ^, 
v T (fl) = qPjy/m^ + (qPn) 2 (n) = MA) 

are used to derive q, m/i and mj^. Another velocity v T (f2) = qP n / \J m j 2 + (?-Pn/c) 2 (/2) = 
Vm{f2) is used to suppress the contribution from the first case. 

To search for events containing a pair of quarks fl and /2 with same sign charge qe and 
(1 — q)e and different masses m/i and m/2 such as ud and us, the events are required to 
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have odd number of charged tracks. The total momentum, energy and one of the velocities 
P r (tot) = qP n (fl) + (1 - g)P n (/2) + others = 0, 
£ r (tot) = y / (gP n ) 2 (/l)+m2 1 + s J{{l-q)P n f{f2)+m} 2 + others = ^ 

v T (fl) = qPjy/m^ + (qP n ) 2 (fl) = v m (fl) 

are used to derive q, m/i and mj 2 . Another velocity v r (f2) = (1 — 
g)P n / ■\Jrn^ 2 + ((1 ~~ Q')-Pn) 2 (/2) = v m (f2) is used to suppress the contribution from the first 
case. 

To search for a pair of quarks / and / with opposite charge qe in events with only two 
charge tracks, 

£ r (tot) = yJ(qP n )*(f)+m} + y/(qP n )*(f)+m 2 f = V~s, 

m f = qP n (fWl-vl(f)/v m (f) 

are used to derive q and m/. Another velocity v T (f) = qP n / \j m2 + (qPn/c) 2 (f) = v m (f) is 
used to suppress the contribution from the first case. 

Search for events containing three free quarks or one free quark and one diquark can 
be done in the similar way. The two simplest examples are A ++ uuu + others and 
/Srddd + others and their charge conjugation, where A decays to A^. 

To estimate the power of the technique above, it can be used to select fractional charge 
particles from a Monte Carlo sample set, which includes a signal sample and a background 
sample. In this article, 200 events of if)(2S) —> ddJ/ip, J/ip — > e + e~ and 200 events of 
ip(2S) uuJ/ip, J/ip — > e + e~ through BESIII detector simulation' 11 ' construct the signal 
sample. The background sample is 40 million events of ip(2S) inclusive decays, in which 
each decay branching ratio is from the particle data group (PDG). d (d) and u (u) carry 1/3 
and 2/3 of the electron charge respectively. Their masses are assumed to be 0.005, 0.025, 
0.045, 0.065 and 0.085 GeV . 

To pick up the signal events and suppress the background events, the events in the sample 
set are selected if they have four charge tracks with zero total charge. The scalar sum and 
the vector sum of the four tracks' momenta are required to larger than 3.8 GeV and larger 
than 0.12 GeV to remove the background events with normal four tracks. Then a pair of 
opposite charge tracks' momenta are adjusted by multiplying a factor of q so that the event 
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satisfy the momentum conservation. The events remain if the adjusted total momentum is 
less than 0.017 GeV, i.e. 

P r (tot) = / Mf) + mG)+Me+)+ Pl (e-)y < 0.017 GeV, 

y i=x, y, z 

Where / and / are supposed to be a pair of opposite fractional charge particles, 
p(/)> P(f)i P( e+ ) an( i P( e ) are measured momenta for the four tracks. If more than 
one pair of opposite charge tracks satisfy the requirement, the combination with minimal 
adjusted total momentum is chosen. The unadjusted pair is supposed to be a pair of e + e~ 
and is required to have invariant mass | m e + e - —3.097 |< 0.100 GeV. The electron (positron) 
track is required to have zero hit in the muon detector. 

Instead to calculate the / and / masses, the masses squared are derived to be m 2 = 
?V(/)(1 - v 2 (f))/v 2 (f), mj2 = gV(7)(l - v 2 (f))/v 2 (f), where v(f) and v(f) are the 
measured velocities and often larger than 1 due to measurement errors. 

Then the total energy is required to be | E'(tot) — 3.686 |< 0.15 GeV, where 

E(tot) = ^(qp) 2 (f)+m 2 + ^(qp) 2 (f) + mL + sj {p) 2 {e+) + m 2 e+ + yfi^P) + m\_ . 

The number of events for the remaining events as well as the measured charge q and 
mass to/ (toj) for the signal channels are listed in table 1 and table 2. The efficiency to 
select the first channel is larger than that of the second channel, because the d and d s 
bending radiuses in the detector are two time those of the u and u if they carry the same 
real momentum and the d and d arrive at the subdetector of the flight time easier. The 
Fig. 1 shows the distribution of the charge q d and Fig. 2 shows the scatter of the mass m d 
{mq) versus the charge qd- It can be seen from the two plots that the measured charge qa 
concentrates at 0.333 and mass m d concentrates 0.010 GeV for the channel ^>(2S) — > ddJ 'ftp 
with the input mass m d = 0.005 GeV. 22 events remain from the background sample, 
which are also plotted in Fig.l and Fig. 2. There is only one event in the d charge window 
| q d — 0.333 |< 0.040 and 4 events in the u charge window | q u — 0.666 |< 0.040. 

In summary, an inclusive sample can be divided into exclusive channels 
riiTvs + niKs + n^ps + n^e+e^s + ns/i + /i~s + n^s + n 7 fs, where / denotes 
a fractional charge particle. For a fractional charge track, its nominal momentum, derived 
by assuming its charge equal to 1 when it is reconstructed, increases by a factor 1/q com- 
pared with its real momentum. Then the nominal total momentum and energy for the event 
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Table 1 The measured charge qa, mass (m^) (GeV) and events for the channel 

^(25) -> ddJji) 



input mass 


0.005 


0.025 


0.045 


0.065 


0.085 


output mass m 2 d 


0.010 2 ± 0.026 2 


0.022 2 ± 0.026 2 


0.037 2 ± 0.031 2 


0.053 2 ± 0.026 2 


0.078 2 ± 0.029 2 


charge q d 


0.339 ± 0.003 


0.333 ± 0.002 


0.333 ± 0.003 


0.329 ± 0.003 


0.323 ± 0.002 


events 


36 


45 


38 


32 


38 



Table 2 The measured charge q u , mass m u {rriu) {GeV) and events for the channel 

ij>(2S) -> uuJ/tP 



input mass m u 


0.005 


0.025 


0.045 


0.065 


0.085 


output mass m 2 


0.016 2 ± 0.031 2 


0.033 2 ± 0.029 2 


0.044 2 ± 0.031 2 


0.064 2 ± 0.029 2 


0.072 2 ± 0.024 2 


charge q u 


0.635 ± 0.006 


0.662 ± 0.005 


0.660 ± 0.004 


0.654 ± 0.009 


0.659 ± 0.004 


events 


28 


24 


33 


32 


27 



containing this kind of tracks will not satisfy the conservation law, which can help find frac- 
tional charge particles described by the two parameters (q, rrif) by analyzing each exclusive 
channel. Between the two parameters, q is more sensitive than m/. Whether the distribution 
of mj concentrates at one or more points very well depends on the momentum resolution 
for charged tracks and the energy resolution for neutral tracks. The effect of | 1 jq — 1 | on 
the momentum for charged tracks is much larger than the resolutions. Let q = 2/3, then 
| l/q — 1 | is 50%, while the momentum resolution is 0.5% for charged tracks and 2.5% for 
neutral tracks for BESIIl' 11 '. The larger | l/q — 1 | is, the more sensitive the method is to 
search for the fractional charge particles, especially for the type of //. Another advantage of 
the method is independence of any theory models. Even if fractional charge particles carry 
continuum masses, the method can be used to search for them, while the dE/dx ionization 
energy loss method is not adequate. 
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FIG. 1: The distribution of the measured charge qd for ip(2S) — > ddJ/ip, J/tp — > e + e~ with the 
input rrid = 0.005 GeV. The shaded area is the signal events. The others are background events. 
The statistics is only for the signal events. 
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FIG. 2: The scatter of the measured m d versus the charge qd for ip(2S) — > ddJ/ip, J/ip 



with the input = 0.005 GeV. The squares are the signal events. The others are background 
events. The statistics is only for the signal events. 
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